Unitarity Corrections to the Structure Functions through the Dipole Picture 
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We study the dipole picture for the description of deep inelastic scattering, focusing on the 
structure functions which are driven directly by the gluon distribution. One performs estimates using 
the effective dipole cross section given by the Glauber- Mueller approach in QCD, which encodes the 
corrections due to the unitarity effects associated with the saturation phenomenon. We also address 
issues about frame invariance of the calculations when analysing the observables. 
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I. INTRODUCTION 



The scattering experiments on deep inelastic electron-proton (DIS) at HERA have provided measurements of the 
inclusive structure function F2{x, Q^) as well as of the Fl and in very small Bjorken variable x values {x ^ 10^^). 
' In these collisions the proton target is analysed through a hard probe with virtuality = —q^. The momentum 
fraction is a; ~ /2p.q, where p and q are the four-momenta of the incoming proton and of the virtual photon 
probe. In the kinematical region of small x, the gluon is the leading parton driving the behavior of the deep inelastic 
^ ' observables. The standard QCD evolution gives a powerlike growth for the gluon distribution and related quantities, 
' and this feature leads, in principle, to the unitarity violation at asymptotic energies, requiring a control of the gluon 
^\ : distribution in high energies. In the partonic language, in the infinite momentum frame, the small momentum fraction 
' region corresponds to the high parton density domain, which is connected with the black disk limit of the proton 
, target and with the parton recombination phenomenon. These issues can be addressed through a non-linear dynamics 
• beyond the usual DGLAP formalism (for a review, see |^). The complete knowledge about the non-linear dynamical 
\^ ' regime plays an important role in the theoretical description of the reactions in the forthcoming experiments RHIC 
^ ' and LHC. 

In the Breit frame, the QCD factorization theorem allows to calculate the scattering processes through the convo- 
lution of the partonic subprocess with the parton distribution functions (pdf's). There, the degrees of freedom of the 
O ' theory are the quasi- free partons (quarks and giuons). On the other hand, more recently the theoretical description 
D of the small x physics has been widely analyzed in the target rest frame, which is a powerful tool concerning an 
unified picture for both the inclusive and the diffractive scatterings, including vector meson production Now, the 
degrees of freedom are the color dipoles, which are the most simple configurations considering the virtual photon Fock 
states expansion. Its main appeal is a quite simplified picture for the different mentioned processes, based on general 
properties of quantum mechanics. 

The description of DIS in the color dipole picture is quite intuitive, allowing a simple representation instead of the 
involved one from the Breit (infinite momentum) frame, and such framework was proposed by Gribov many years 
ago Q. Considering small values of the Bjorken variable x, the virtual photon fluctuates into a qq pair (dipole) with 
fixed transverse separation r at large distances upstream of the target and interacts in a short time with the proton. 
More complicated configurations should be considered for larger transverse size systems, for instance the photon 
Fock state qq + gluon. An immediate consequence of the lifetime of the pair {Ic = l/2mpx) to be bigger than the 
interaction's one is the factorization between the photon wavefunction and the cross section dipole-proton in the 7* p 
total cross section. The wavefunctions are perturbatively calculable, namely through QED for the qq configuration |5[| 
and from QCD for qqG The effective dipole cross section should be modeled and encloses both perturbative and 
non-perturbative content. However, since the interaction strength relies only on the configuration of the interacting 
system the dipole cross section turns out to be universal and may be employed in a wide variety of small x processes 
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Currently, there are several models for the dipole cross section based either in pure phenomenological parametriza- 
tions or in a more theoretical ground (for a review, see Ref. |^). The main feature in those models is the description 
of the energy dependence of the interaction, namely taking into account the interplay between hard and soft domains. 
The dipole cross section should be consistent with the sharp growth on energy at small transverse separation r (large 
gluon density) and a softer behavior for larger r (Regge-like phenomenology). An additional ingredient is the expec- 
tation for saturation effects in high energies as a consequence of the unitarity requirements. Indeed, the growth of the 
gluon distribution should be tamed at very small x and it has been found that the corrections are important already 
in the present HERA kinematics 0. Moreover, such effects are associated with higher twist contributions concerning 
the standard linear evolution equation, i.e. the DGLAP formalism 

Here, we take into account a formalism providing the unitarity corrections to the deep inelastic scattering at small 
X, namely the Glauber-Mueller approach in QCD. It was introduced by A. Mueller who developed the Glauber 
formalism to study saturation effects in the quark and gluon distributions in the nucleus considering the heavy onium 
scattering. Afterwards, the authors of Ref. extended that approach that has as a limit the GLR results ||ll|. It is 
obtained an evolution equation taking into account the unitarity corrections (perturbative shadowing), generating a 
non-linear dynamics which is related with higher twist contributions. Its main characteristic is to provide a theoretical 
framework for the saturation effects, relying on the multiscattering of the pQCD Pomeron. The latter is represented 
through the usual gluonic ladder in the double logarithmic approximation. 

Summarizing the Glauber-Mueller approach, the gluonic content of the nucleus or nucleons is obtained in the 
following way: in the rest frame, a virtual probe (gluon) decays into a gluon pair interacting with the nucleon inside 
the nucleus. The multiple scatterings of the pair give rise to the unitarization of the corresponding cross section. 
The calculations are performed in the double logarithmic approximation (DLA), corresponding to the condition 
Q^s ^ 7g Ij where and ^'^e the QCD coupling constant and the gluon anomalous dimension, respectively. 
In this approximation, the transverse separation of the pair remains fixed allowing an eikonal description for the 
interaction gluon pair- nucleon through incoherent multiscatterings pO[ | . The cross section for the interaction can be 
expressed in terms of the nucleon gluon distribution xG{x, Q"^) and of the transverse separation r of the gluon pair. 
The procedure for an initial state quark-antiquark pair is similar, up to proper color coefficients. Such a formulation has 
produced comprehensive phenomenological applications: the inclusive structure function F2 | p2| , the longitudinal 
and charmed ones , the logarithmic slope dF2 / d log and other related quantities have been calculated 

in the Breit system. The respective nuclear case has been investigated in Refs. [|l^ and the equivalence with other high 
density QCD approaches has been reported in Refs. ||l6|] . Moreover, the asymptotic limit of the inclusive structure 
function and its logarithmic slope are estimated in [ p.7| . Regarding the rest frame, Glauber-Mueller has been also 
used to estimate the saturation effects for DIS and diffractive dissociation as well as it is considered as initial 
condition for a high energy evolution equation (for a review, see |^]). Therefore, the Glauber-Mueller approach gives 
a good framework for the unitarity effects (saturation) in the nucleon and nuclear sectors, providing the dynamics of 
the observables in a quantitative level. 

In this work we make use of the parton saturation formalism to study the description of the observables driven by 
the gluonic content of the proton in the color dipole picture. The inclusive structure function F2 is calculated properly, 
disregarding approximations commonly considered in previous calculations |T2|jr^ , [l7t . The structure functions Fl and 
^2'^ are presented for the first time using the Glauber-Mueller approach and the rest frame in comparison with the 
experimental data. The saturation effects are included in the effective dipole cross section corresponding to the small 
size dipole contributions. The large dipole sizes are taken into account through an ansatz for the non-perturbative 
region. Here, we choose to freeze the gluon distribution at large distances jist . The dipole framework provides a clear 
identification in the transverse distance r range where the perturbative and non-perturbative sectors (soft domain) 
contribute to the calculated quantities. Large dipoles correspond to soft domain and the small ones are connected 
with the hard sector. As we will see, the photon wavefunctions play the role of a weight function selecting small 
dipole sizes, with a non-negligible contribution coming from the large dipole sizes. However, when considering the 
production of heavy quark pairs (charm and bottom), their masses largely diminish the soft contribution. These issues 
are addressed throughout this paper. 

This work is organized as follows. In the next section we shortly review the deep inelastic scattering in the rest 
frame, introducing the main formulae for the further analysis. In Sec. (3) one addresses the effective dipole cross 
section considering saturation effects (unitarity corrections) encoded in the Glauber-Mueller approach, pointing out 
its main properties and discussing the large transverse separation contribution. The Sec. (4) is devoted to calculate 
the theoretical estimates for the HERA observables, focusing on those ones dominated by the gluon content. In the 
last section we draw our conclusions and comments. 
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II. THE DEEP INELASTIC SCATTERING IN THE PROTON REST FRAME 



In the reference frame where the target (proton) has infinite momentum the usual description of the dynamics 
evolution is the following: the emitted partons from the target remain quasi-free for enough time in such a way that 
the virtual probe (photon) detects them as real particles (asymptotic freedom). On the other hand, in the small x 
region it is suitable to describe the evolution in a system where the target is at rest: in this situation the evolution 
is related to the partonic fluctuations of the probe and their interactions with the nuclcon. The rest frame physical 
picture is advantageous since the lifetimes of the photon fluctuation and of the interaction process are well defined [p^ . 
The simplest case is the quark-antiquark state (color dipole), which is the leading configuration for small transverse 
size systems. Its life time can be estimated by the uncertainty principle through the energy fluctuation associated with 
the emerging pair. The well known coherence length is expressed as Ic = l/{2xmp), where x is the Bjorken variable 
and rup the proton mass. For instance, in deep inelastic at HERA kinematics reaching at x ^ 10~^, the coherence 
length is about 10^ fm, which is a distance larger than the radius of any atomic nuclei. An immediate consequence 
of such a picture is the factorization between virtual photon wavefunctions and the interaction cross section on the 
corresponding amplitude in the impact parameter space representation. 




FIG. 1. Representation of the deep inelastic process, where the upper blob corresponds to the photon impact factor and the 
bottom one represents the proton impact factor. The corresponding two first orders in perturbative expansion (BFKL-like) are 
depicted. 

A striking consequence of the formulation above is that the photoabsortion cross section can be derived from the 
expectation value of the interaction cross section for the multiparticle Fock states of the virtual photon weighted 
by the light-cone wave functions of these states ||]. The scattering matrix has an exact diagonalization in the 7* 
Fock states representation, where the partial-wave amplitudes subjected to the s-channel unitarization are identified 
properly. Then, the photoabsortion cross section can be cast into the quantum mechanical factorizcd form, 

allix, Q2) ^ J d\ dz \^tA^^ i-)|2 a'i'P°i-(x,z,r) , (1) 

The formulation above is valid even beyond perturbation theory, since it is determined from the space-time structure 
of the process. The 5't,l(-z, r) are the photon wavefunctions (for transverse, T, and longitudinal, L, polarizations) 
describing the pair configuration, where z and 1 — z are the fraction of the photon's light-cone momentum carried by 
the quark and antiquark of the pair, respectively. The transverse separation of the pair is r. The precise normalizations 
of the wavefunctions can be determined through the Fock expansion |7phys >= |7baro > + ■N'qq |'?^>, where Z3 is 
the 7-wavefunction renormalization constant, |7baic > denotes the bare state and Afqq is the coefficient determining 
the probability of the qq pair fluctuation in the photon [^0| . Considering completely normalized states, then A/'gg — 
1 — Z3 and the normalization for transverse photons is obtained from J dz(fr |^'t.l(z, r)p — M^q. The remaining 
normalizations, i.e., longitudinal component and cross section, are consequently fixed [ pO| . The explicit expressions 
are well known, 

nf 

(z, r) p = ^ ^ { [ + (1 - zf ] Kf (e r) + KI {e r) } (2) 

i 
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and 



i 

Clarifying the notation, the auxihary variable = z(l — z)Q^ + m^, with the light quark mass, and Kq and 
Ki are the Mc Donald functions of rank zero and one, respectively. 

The quantity (t'^'p°'''(x, z, r) is interpreted as the cross section of the scattering of the effective dipole with fixed 
transverse separation r H . It is directly dependent on the unintegrated gluon distribution, 

J J ki z 

The unintegrated gluon distribution, J-{ ^, ki), vanishes at the gluon transverse momentum |ki| (and similarly 
for the factor (1 — e'*'^''")) due to gauge invariance, and therefore the dipole cross section has to be infrared finite. The 
most important feature of the dipole cross section is its universal character, namely it depends only on the transverse 
separation r of the color dipole. The dependence on the external probe, i.e. the photon virtuality, is included in the 
wavefunctions. 

The main technical difhculty in Eq. is to model the unintegrated gluon distribution function in a suitable way, 
mainly in the small transverse momentum (fc^) region (infrared sector). In particular, to obtain these distributions 
one should solve numerically an evolution equation, which makes the procedure cumbersome for practical use (instead, 
for a prompt parametrization see ||2l|] ). In general, this is avoided introducing an ansatz for the effective dipole cross 
section and analyzing the process in the impact parameter space. The main feature of the current models in the 
literature is to interpolate the physical regions of small transverse separations (QCD-parton improved model picture) 
and the large ones (Regge-soft picture) . Below, we quote two of them, which have connections with the work performed 
here. 

The phenomenological saturation model of Golec-Biernat and Wiisthoff (hereafter GBW) gives a good description 
of DIS data |23j. The corresponding cross section interpolates between color transparency, i.e. (j'^'p°'° ^ at small 
r, and constant cross section ctq at large r (confinement). Such a procedure ensures Q^-saturation, while parton 
saturation at low x is obtained with an eikonal-inspired shape for the dipole cross section, 

a«W(x,r^) = ao(l-e-^^/«-(^)) , (5) 
Wa \ ■''0 

where the parameters cto = 23.03 mb, xq = 3.04 10"'' and A = 0.288 are fitted to the HERA DIS inclusive data 
with X < 10~2, whereas Qq = 1 GeV^ sets the dimension. The x-dependent saturation radius Rq{x) scales the 
pair separation r in the cross section and is associated with the mean separation between partons in the nucleon. 
This approach was used to describe diffractive dissociation in a parameter-free way, considering also the required 
qqG configuration. Some criticism to the GBW model, mainly concerning a better knowledge of the parton (gluon) 
distribution at low r and its saturation at small x are postponed to the next section. The most of them can also to 
be found in Refs. 

Our work in this paper is closer to the McDermott et al. one [Q (hereafter McDFGS), where perturbative QCD 
relates the dipole cross section to the leading logarithmic gluon distribution function in the proton, at LLA{Q^) 
accuracy, 

^McDFGS ^!l!^fiQ!) ^2 ^G(x,Q2), (6) 
o 

where the identification w is allowed at leading-log level. A study of the dipole cross section is carried for 
all transverse separations r in p^ . The small r region is described by the usual gluon pdfs, evoluting with the 
scale = lO/r^, while for large separations (rpion > 0.65 fm) the dipole cross section is driven by the pion-proton 
contribution with the typical soft energy behavior from the hadronic sector. Moreover, the taming of the parton 
density is implemented by hand starting at the named critical transverse separation rcrit, stated when the dipole 
cross section reaches one half of its maximum value labeled by the pion-proton cross section. In connection with the 
present work, the Glauber-Mueller approach gives the Eq. (||) at the Born level. Instead of using an ad hoc control 
of the gluon distribution, the Glauber-Mueller provides corrections required by unitarity in an eikonal expansion. For 
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the large r region, we choose to follow a similar procedure from the GBW model, namely saturating the dipole cross 
section (r- independent constant value). 

Having defined the notation and reviewed the main properties settled by the rest frame representation of deep 
inelastic process, in the next section we address the unitarity corrections formalism contained in the Glauber-Mueller 
approach, that we will consider in the calculation of F2, Fl and i^2^^. 



III. THE GLAUBER-MUELLER APPROACH 



The Glauber formalism concerns mainly interactions with a target nucleus, allowing to calculate the amount of 
unitarity corrections to the nuclear cross section. However this approach can be extended to take into account 
the evolution of the partonic densities (saturation) through the multiple scatterings. Below, we review shortly the 
main properties of the Glauber formalism in QCD, either in the nuclear case as its application to the nucleon case. 



We indicate the original papers |10| for a complete presentation. It should be stressed that saturation effects mean 
unitarity corrections to the observables. Indeed, the asymptotic calculations have produced an unified ln(l/x) pattern 
for the cross section and gluon function instead of a effectively saturated one [|7| . 

Since the small x limit is driven by the gluonic interactions, we consider a virtual probe G* with invariant mass 
which decays into a gluon pair GG having a transverse separation r and transverse momentum k. Then, the pair 
interacts with the target through a gluonic ladder at fixed transverse separation. Following the discussion from the 
previous section, the photoabsortion cross section for the probe particle in terms of x and virtuality in the nuclear 
case is 

Q2) ^ j J^' dz \^gg{z, r,Q2)|2 a«G-„ucieus(^^ ^.^ ^ 

where the variables have the same identification as in the qq pair discussion. The quantity ^'gg is the light-cone 
wavefunction for the gluon pair. The Glauber's multiscattering theory employs the phase shift method to describe 
processes at high energy for an incident particle undergoing into successive scatterings. When x is small, the coherence 
length Ic is bigger than the mean radius Ra for all nuclei and the interaction of the initial parton stands through the 
entire nuclear path, providing coherent interactions with all target partons along the distance Ic- These scatterings 
are coherent in an interaction of a hadronic fluctuaction, i.e. quark or gluon pairs, with the nucleus. There are 
interference effects among them, generating a reduction in the nuclear cross section, cr^ucicus ^ ^^nucicon^ Otherwise, 
when Ic < full incoherent scatterings occur leading to the expectation that the nuclear cross section equals to 
^ ^nucleoli _ rpj^g ^gjj known Glauber formula for the total cross section of a hadronic state with the nucleus is 

^nucleus ^2 [ d^h (l~ '^"-loon SA(b)A ^ (8) 



where SaO^) is a profile function containing the dependence on the impact parameter b, which is the conjugate 
variable to the momentum transfer t. It is related to the nucleon distribution inside the nucleus and encodes the 
information about the angular distribution of the scattering. We discuss its particular shape later on. The Eq. (|^) is 
quite general and allows to describe the hadronic fluctuations of the photon or any virtual probe, as saw in the last 
section. 

Considering the scattering amplitude dependent on the usual Mandelstan variables s and t, now written in the 
impact parameter representation b, 

a(s,b)^i- y'd2qe-''ib^(s,t = -q2). (9) 

the corresponding total and elastic cross sections (from Optical theorem) are rewritten in the impact parameter 
representation (b) as 

atot = 4:TrImAis,0) = 2 J d^hlm a{s,h) , (10) 

a,, = J d^h\a{s,h)\\ (11) 

A very important property when treating the scattering in the impact parameter space is the simple definition for 
the unitarity constraint UM- At fixed b, the constraint can be expressed in the following way, 
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CTtot = (^el + (^inel , (12) 

2 Im a{s, b) = \a{s, b) ^ + b) , (13) 

with Cm(s,b) denoting the sum of contributions from all the inelastic channels. The constraint above has a simple 
solution. If the real part of the scattering amplitude vanishes at the high energy limit, corresponding to small x 
values, the solution is 

a(s,b) h -e-5f^ 1 , (14) 



atot = 2 I d^h 



1 _ g-^n(s,b) 



(15) 



where the opacity f2 is an arbitrary real function and it should be determined by a detailed model for the interaction. 
The opacity function has a simple physical interpretation, namely e~^ corresponds to the probability that no inelastic 
scatterings with the target occur. To provide the connection with the Glauber formalism, the opacity function can 
be written in the factorized form il(s,b) = il{s) S{h), considering S'(b) normalized as / d^bS'(b) = 1 (for a detailed 
discussion, see i.e. p5[|). 

From Eqs. (jl^) and (^), we identify the opacity Q{s « Q'^/x;r) — cr""'^'°°"(a::, r). Moreover, it has been found that 
the same formalism for multiple scatterings can be applied to the nucleon case. To proceed, we should determine the 
GG cross section. The gluon pair cross section is equivalent to the quark pair one, up to a color factor {a'^'~^ — |o'''^). 
The {qq pair) dipole-proton cross section is well known |l^j2^, calculated starting from Eq. (^), and in double 
logarithmic approximation (DLA) has the following form 

-Zlcon(-,0 = ^^^^^-G(.,Q^) (16) 



with the r-dependent scale Q — r^/r . Considering Eq. ( |16D one can connect directly the dipole picture with the 
usual parton distributions (gluon), since they are solutions of the DGLAP equations. In our case, we follow the 
calculations in Ref. (l^,^^ and consider the effective scale = Such a value differs from |^^, where it is 

r§ — 10, which is obtained by an averaging procedure on the transverse size integral of -F^. However, in further 
studies in vector mesons it was found that Tq ranges from 4-15, and F2 and are not sensitive to those variations. 
Thus, these values are consistent in leading logarithmic approximation. 

From the above expression, one obtains a dipole cross section satisfying the unitarity constraint and a framework 
to study the unitarity effects (saturation) in the gluon DGLAP distribution function. Hence, hereafter we use the 
Glauber-Mueller dipole cross section given by 



GM 
^dipole 



j cfh (l-e-^-lVo„('^.r)s(b)^ (^7) 



Here, some comments are in order. The Glauber-Mueller approach is valid in the small x region, and the gluon 
emission is described in the double logarithmic approximation of perturbative QCD. The interaction of the quark 
pair with the nucleon (proton) occurs through ladder diagrams exchange, satisfying the DGLAP equations in the 
DLA limit. The high energy limit allows to treat successive scatterings as independent collisions, meaning the process 
described by the eikonal picture of a relativistic particle crossing the nucleus. Moreover, as a consequence of no 
correlation among nucleons inside the nucleus (in the nuclear case), there are no correlations among partons from 
different partonic cascades, stressing that only the fastest parton interacts with the target. The corrections coming 
from the slowest partons in the cascade (emitted by the pair) lead to the AGL non- linear evolution equation [ p^ , and 
they have been considered recently to describe diffractive DIS in Ref. |lq| . Regarding criticisms to the Glauber-Mueller 



approach, we indicate the recent paper |26| for a complete discussion about the eikonal- like models, concerning their 
advantages and limitations as well as pointing out the improvements to be taken into account to introduce the proper 
corrections. 

Now, we proceed to calculate numerical estimates of the dipole cross section using the Glauber-Mueller approach 
through Eqs. Then we are calculating saturation effects in the color dipole picture. Firstly, we need to 

discuss the profile function S{b). This function contains information about the angular distribution in the scattering, 
namely the t dependence (quark pair-ladder and proton-ladder couplings) . Both of them can be approximated by an 



exponential parametrization, leading to a simple gaussian shape in the impact parameter space, S{h) = e 

where A is the atomic number and Ra is the target radius. We will keep this notation although we are only concerned 
with the nucleon case. The R\ value should be determined from data, ranging between 5—10 GeV^^ for the proton 



A 
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case (see discussions in [|18[|7|]). For nuclear reactions, a more suitable shape for the profile should be taken into 
account, since the gaussian approximation is no longer appropriate to describe the nuclear profile for large A. Here, 
we have used the value ( R\ = 5 GeV~^) obtained from a good description of both inclusive structure function and 
its derivative jl^. Such a value corresponds to significative unitarity corrections to the standard DGLAP input even 
in the current HERA kinematics. 

Now, we discuss in a detailed way the main characteristics emerging from the dipole cross section Eq. (|l7|). In 
order to do this, in Fig. (||) one shows the Glauber-Mueller dipole cross section as a function of dipole transverse 
size r at fixed momentum fraction x. For a better illustration on the expected partonic saturation effects, we run the 
Bjorken x down to a quite small value x = 10~^ (THERA region). Hereafter, we are using the GRV gluon distribution 
at leading order |27| in the input Eq. (^6|), whose choice we justify below. The solid lines correspond to the dipole 
cross section calculation, Eq. (p^, whereas the dashed lines are the GBW model presented for comparison. The 
general shape in terms of the dipole size comes from the balancing between the color transparency crdip ~ r'^ behavior 
and the gluon distribution. 

Here some comments about the large transverse separation are in order: although perturbative QCD provides 
reliable results at small distances (small dipole sizes), the nonperturbative sector is not still completely understood. 
The usual pdf's are evoluted from a perturbative initial scale Qq ^ 1 GeV^, and one has little information about 
the behavior at < Qq. In general one makes use of Regge phenomenology to estimate those contributions (see, 
for instance p4|). Thus, extrapolating to lower virtuality regions (large dipole sizes) one needs a parametrization 
regarding the nonperturbative sector. 

This is the main justification of the use the GRV94 parametrization ||2^ in our calculations. Bearing in mind that 
= its evolution initial scale is Qq = 0.4 GeV^ allowing to scan dipole sizes up to Vcut = ^ GeV~^ (= 

0.62 fm). For the most recent parametrizations, where Q\ ^ \ GeV'^ (r^ut ~ 0.4 fm) the amount of nonperturbative 
contribution in the calculations should increase. An additional advantage is that GRV94 does not include non-linear 
effects to the DGLAP evolution since the parametrization was obtained from rather large x values. This feature 
ensures that the parametrization does not include unitarity corrections (perturbative shadowing effects) in the initial 
scale. 

Now, we should introduce an ansatz for the large transverse separation region. A more phenomenological way is 
to match the pQCD dipole cross section with the typical hadronic one fT,rAr at rcut as performed in |^^. However, 
due to the significant growth of the pQCD dipole cross section at high energies, we choose an alternative ansatz: the 
gluon distribution is frozen at scale rcuti namely xG{x, Qcut)- Then, for the large distance contribution r < r-cut the 
gluon distribution reads as 

xG{x,Q^ <Ql)^%xGix,Q^ ^Ql), (18) 

leading to the correct behavior x G ^ as ^ 0. In a more sophisticated case, one can substitute the freezed 
scale Qcut by the saturation scale Qg{x) to take into account a realistic value of the gluon anomalous dimension in 
all kinematic region (see correlated issues in p^). 

Recently, the phenomenological model of Ref. |Q has produced a good description of HERA data in both inclusive 
and diffractive processes. It is constructed interplaying the color transparency behavior adip ~ r'^ at small dipole 
sizes and a flat (saturated) behavior at large dipole sizes adip ctq (confinement). The expression has the eikonal-like 
form, 



1-cxpf ■ ^0 



4{x/xq)^ 



(19) 



where Qq = 1 GeV^ and the three fitted parameters are ctq = 23.03 mb, xq = 3.04 10^** and A = 0.288 and the notation 
for the saturation radius Ro{x) = {x/ Xq)'^^''^ . The GBW total cross section lies below the typical hadronic cross section: 
for instance in the pion-proton case, convoluting the pion wavefunction squared with the GBW dipole cross section 
we would have a constant cross section at high energies (jfJt^ — '^tot' ■ Despite describing data in good agreement, 
GBW has some details that deserve some discussions: the approach does not present a dynamical hypothesis for 
the saturation phenomena and does not match DGLAP evolution. In GBW, saturation is characterized by the x- 
dependent saturation radius Ql{x) = 1/Rq{x) instead of the scale coming from Glauber-Mueller, kc{x,QI) = 1, 
which can be easily extended for the nuclear case pO|. 
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FIG. 2. The dipole cross section as a function of r at fixed x 
GBW model results (dashed curves) are also included. 
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The Glauber-Mueller corresponds to the solid line and the 



In Fig. (U) one shows the Glauber-Mueller dipole cross section and the GBW model [g3| . We choose to compare them 
due to the fact that GBW is actually a particular case of the Glauber-Mueller approach, considering an oversimplified 
profile function. We should show below that this fact allows to construct an extended saturation model with DGLAP 
evolution (BGK model). In the lower plots, where x = 10~^ and x = 10~'^, the GM cross section underestimates 
the GBW one. However, as x decreases the gluon distribution in the proton rises producing a bigger dipole cross 
section. This feature is clear in the upper plots, for smaller x = 10~^ and x = 10~^, where GM overestimates GBW 
by a significant factor mostly at intermediate r values. 

Finally, we show the connection between the Glauber-Mueller approach with the saturation model with DGLAP 
evolution |28| , observing that it is timely since the latter is quite efficient in describing data and its qualitative success 
corroborates quantitative QCD studies. Considering the particular case of central collisions, namely scattering at 
impact parameter b = {S = A/tt R\), the Glauber-Mueller approach produces. 



GM 
dip 



' (t 



0) 



d^h 1 



(20) 



The integration over b can be promptly carryied out, and introducing the notation for the proper normalization for 

the dipole cross section, ctq = 2 vr J^ '^ db^ = 2tt R\ , the Eq. (|20| ) recovers the simple expression for the saturation 
model DGLAP evoluted |28|, 



c^dipofo(a^>i") = cTo 



1 — exp — 



TT^ asifJ-^) xG{x, fi'^) 

3 CTo 



(21) 



For a phenomenological analysis, the parameter ctq and the scale are determined from data in [ p8| . In our case, it 
assumes the well defined value CTq = 12.22 mb, using R\ = 5 GeV~^. For a larger radius, for instance R\ = 10 GeV~^, 
one obtains a value ctq = 24 mb, closer to the GBW one. Here, the virtuality scale is = 4/r^), whereas BGK 
choose the parametric form = fi^ = C/r^ + /ip and a two-parameter initial condition for the gluon distribution 
function. 



IV. UNITARITY EFFECTS IN EP COLLISIONS 



This section is devoted to the study and estimate of the gluon driven observables measured at HERA kinematical 
domain in the rest frame. The first one is the inclusive structure function F2(x, Q^), the main quantity testing the small 
X physics. The unitarity corrections are well established for this observable considering Glauber-Mueller approach |l^ ] 
as well as its high energy asymptotics namely the black disk limit. We review these issues considering the dipole 
picture (rest frame), using a more complete analysis similar to |p3| , p4| , but mostly, discuss in detail the role played 
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by the nonperturbative physics needed to describe the structure function, and where in the transverse separation r 
range it starts to be important. 

The longitudinal structure function Ft,{x , Q^) is also addressed, verifying the frame invariance in comparison with 
previous laboratory frame calculations [|l3| . The longitudinal wavefunction strongly suppress large r contributions, 
thus selecting smaller nonperturbative contribution in comparison with the F2 case. Moreover, is one of the 
main observables scanning possible higher twist corrections in the standard Operator Product Expansion (OPE) 
1^. Therefore, a reasonable description of this quantity suggests that the Glauber-Mueller formalism (or similar 
eikonal-like approaches) take into account the most important contributions to the complete higher-twist corrections 
at current kinematical regimes. 

The structure function i^2'^(a;, Q^) gives the charm quark content on the proton and is directly driven by the gluon 
distribution. Therefore it is a powerful observable to scan saturation effects in the small x region. However the current 
experimental status requires more dedicated measurements and a better statistics. We verify a consistent description 
in the rest frame corroborating the similar analysis in the dipole models p^j29[ | and in those ones considering unitarity 
corrections in the laboratory frame |0| . 



A. The Inclusive Structure Function F2{x,Q^) 

Now we perform estimates for the inclusive structure function in the rest frame considering the Glauber-Mueller 
dipole cross section. The expression for F2, with the explicit integration limits on photon momentum fraction z and 
transverse separation r is. 



ITT aom Jo J a 



(22) 




FIG. 3. The Glauber- Mueller (GM) result for the F2(x,Q^) structure function. It is shown the transverse contribution 
(dot-dashed), the longitudinal one (dashed) and total one (solid line). One considers light quarks, target radius R\ = 5 GeV~^ 
and frozen gluon distribution at large r > rcut. The input gluon distribution is GRV94LO [27] and data from HI Collaboration 
[30]. 



The notation has been introduced in the previous sections. In the Fig. (^) one shows the estimates for the structure 
function for representative virtualities from the latest HI Collaboration measurements [^ . The longitudinal and 
transverse contributions are shown separately, the longitudinal one being subdominant as is well known. An effective 
light quark mass (u, d, s quarks) was taken, with the value niq = 300 MeV, and the target radius is considered R\ = 5 
GeV^^, in agreement with Ref. ||l4|. It should be stressed that this value leads to larger saturation corrections rather 
than using radius ranging over ^ 8 - 15 GeV~2. The soft contribution comes from the freezing of the gluon 
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distribution at large transverse separation as discussed at the previous section. The gluon distribution considered is 
GRV94 at leading order xG'''^{x, whose choice has been justified in the previous section. 

From the plots we verify a good agreement in the normalization, however the slope seems quite steep. This fact 
is due to the modeling for the soft contribution and it suggests that a more suitable nonperturbativc input should 
be taken. Indeed, in Ref. |24[ such a question is addressed, claiming that the correct input is the pion-proton cross 
section parametrized through the Donnachie-Landshoff pomcron. It is found that the large transverse separations 
give a larger contribution at low Q^, whereas it vanishes concerning higher virtualities. 




FIG. 4. The integrand Hggrir, x, Q^) as a function of r for = 2 and = 30 GeV^ at fixed 10" 

To clarify this issue, we can calculate the integrand Hqq[r,x,Q^) used in the r-integration for the cross section 
(Tt = /q°° drHqqT{r, x, Q^), that should have significant nonperturbativc content. We plot in the Fig. (^) this quantity 
at low virtuality = 2 GeV^ and at higher one = 30 GeV^ for momentum fraction ranging on 10^** < a; < 10^^, 
verifying that the main contribution comes from an asymmetric peak at r w 0.15 fm for — 2 GeV^, while it is 
shifted to r « 0.07 at = 30 GeV^. In our calculation, the pcrturbative contribution holds up to rcut — 0.62 
fm, therefore the region r > rcut gives a non marginal contribution to the cross section at low virtualities. Indeed, 
one has found that it reaches about 10 % at = 2 GeV^ and that when the virtuality increases the contribution 
gradually vanishes. In fact, using the most recent pdf's this situation is more critical since rcut is smaller (Qq ~ 1 — 2 
GeV^). This suggests that the photon piece |\E'T(2:,r)p multiplying the dipole cross section enhances the r-integrand 
to smaller r at high Q^, corroborating a similar conclusion already found in Ref. p^ ]. 

To clarify the role played by the soft nonperturbative contribution to the inclusive structure function and to verify 
the frame invariance of the approach, in the Fig. we plot separately the pcrturbative contribution and parametrize 
the soft contribution introducing the nonperturbative structure function = Csoft (1 — ||l^, which is 

added to the pcrturbative one. The soft piece normalization is Csoft = 0.22. Such procedure is done in order to 
compare explicitly the results found in Ref. iQ. Accordingly, we have used just shadowing corrections for the quark 
sector, taking into account only the transverse photon wavefunction and zero quark mass. The integration on the 
transverse separation is taken for < ''^ < l/*3o' with Qq — 0.4 GeV'^ for leading order GRV94 gluon distribution. 
This leads to a residual contribution absorbed in the the soft piece coming from transverse separations r^ < 



We considered the target radius being R\ = 5 GeV^ (supported by |14|), which produces a correction more important 
than the value i?^ = 10 GeV^ It is again verified that the soft contribution is important at small virtualities and 
decreasing as it gets larger. In the plots, the dot-dashed lines represent only the pcrturbative calculations using the 
particular approximations indicated above, and the solid lines represent the results when we add the soft term. 
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FIG. 5. The Glauber-Mueller prediction for the F2 structure function in the rest frame. For sake of comparison with the 
results in [12] , one uses the quark sector (ii^ — 5 GeV~^, niq — 0) and only the transverse wavefunction. Radius integration 
1/Q^ < < 1/Qo s-nd soft Pomeron added (parametrizing the large pair separation - F|°** — Csoft a;~''"*(l — xY^). 

Concluding, we have a theoretical estimate, i.e. no fitting procedure, of the inclusive structure function F2{x,Q^) 
through the Glauber-Mueller approach for the dipole cross section, detecting a non negligible importance of a suitable 
input for the large dipole size region. 




B. The Longitudinal Structure Function Fl{x,Q^) 

As we saw in the Sec. (||), the inclusive structure function can be expressed in terms of the cross sections (Tt and 

ctl for the absortion of virtual photons transversally and longitudinally polarized, ^2(2;, Q^) = — (ctt + o'l)- At 
small X, the longitudinal structure function is 

FLix,Q^) = -^aLix,Q^). (23) 

From QED, the longitudinal photons have zero helicity {h = 0) and therefore they have a virtual character. In 
the naive parton model, the helicity conservation for the electromagnetic vertex implies to the Callan-Gross relation 
F2 = 2xFi and consequently a vanishing value for the longitudinal structure function Fl = F2 — 2xFi , considering the 
scattering photon-quarks (spin 1/2). From QCD theory, this quantity has a non-zero value due to the gluon radiation, 
as is encoded in the Altarelli-Martinelli expression (see discussion in [jisll ) 



27r 



\F2{x,Q^) + ^yG{y,Q^) ( 1 - ^ 
3 9 \ V 



(24) 



where y = Q"^ j sx is the inelasticity variable. Therefore, the structure function Fj^ is an auxiliary observable to detect 
saturation effects (unitarity corrections) in the gluon distribution. 

Experimentally, the determination of the is quite limited, providing few data points. Most recently, the HI 
Collaboration has determined the longitudinal structure function through the reduced double differential cross section 

Gr^F2{x,Q')-^-FL{xM^), (25) 

where F+ = 1 -I- (1 — y)^. For large inelasticity, the reduced cross section becomes (F2 — Fl) and the contribution 
of Fl is enhanced for large y. The longitudinal structure function should be obtained only in the region of large 
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inelasticity, covered in a large range at HERA. In Ref. |3^], two methods were used to perform the extraction: (i) 
for large > 10 GeV^, is obtained through the extrapolation method, using a NLO DGLAP QCD fit (in the 
restrict kinematic range y < 0.35 and < 3.5 GeV^) to extrapolate F2 into the high y region, (ii) At low < 10 
GeV^, the behavior of F2 as a function of \ny is obtained using the derivative method, based on the cross section 
derivative (9(7^/9 In ?;)g2. The data points obtained arc consistent with the previous measurements, however they are 
more precise and lying into a broader kinematical range. Therefore, in the following we use only the new data points 
to analyse. 
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FIG. 6. The Glauber- Mueller estimates for the Fl structure function. One uses light quarks {niq — 300 MeV), target size 
R\ — 5 GeV"'^ and frozen gluon distibution at large r. Data from HI Collaboration [30]. 



In Fig. (g) we present the estimates for the Fl structure function, in representative virtualities as a function of 
X. For the calculations, it was considered light quarks (m, d, s) with effective mass niq — 300 MeV and the target 
radius R\ — 5 GeV~^. The large r region is considered by the freezening of the gluon distribution at that region. 
Our expression for the observable is then. 



47r2 



dipolc V 



(26) 



The behavior is in agreement with the experimental result, either in shape as in normalization. A better description 
can be obtained by fine tunning the target size or the considered gluon distribution function, however it should be 
stressed that the present prediction is parameter-free and determined using the dipole picture taking into account 
unitarity (saturation) effects in the effective dipole cross section. We verify that the rest frame calculation, taking 
into account the dipole degrees of freedom and unitarity effects produces similar conclusions to those ones using the 
Breit system. For instance, in a previous work IT^ , the unitarity corrections to the longitudinal structure function 
were estimated in the laboratory frame considering the Eq. (p^), with unitarized expressions for F2 and xG{x, Q^), 
obtaining that the expected corrections reach to 70 % as ln(l/a;) — 15, namely on the kinematical sector of the 
upcoming THERA project. 

The higher twist corrections to the longitudinal structure function have been pointed out. For instance, Bartels et 
al. 1^ have calculated numerically the twist-four correction obtaining that they are large for Ft and F^, although 
with opposite signs. This fact leads to remaining small effects to the inclusive structure function by almost complete 
cancellation between those contributions. The higher twist content is analyzed considering the model as initial 
condition. 



12 



Concerning F^, it was found that the twist-four correction is large and has negative signal, concluding that a 
leading twist analysis of Fj^ is unreliable for high and not too small x. The results are in agreement with 
the simple parametrization for higher twist (HT) studied by the MRST group in Ref. ||3|], where F^^{x,Q^) = 

F2"'" {x,Q^){l + ^^Qi^^ ). The second term would parametrize the higher twist content. In our case, the unitarity 
corrections provide an important amount of higher twist content, namely it takes into account some of the several 
graphs determining the twist expansion (for recent discussions in these issues, see Q). 



C. The Charm Structure Function Ff{x,Q'^) 

In perturbative QCD, the heavy quark production in electron-proton interaction occurs basicaly through photon- 
gluon fusion, in which the emitted photon interacts with a gluon from the proton generating a quark-antiquark 
pair. Therefore, the heavy quark production allows to determine the gluon distribution and the amount of unitarity 
(saturation) effects for the observable. In particular, charmed mesons have been measured at deep-inelastic at HERA 
and the corresponding structure function i^|'^(a:, Q^) is defined from the differential cross section for the cc pair 
production. 



dx dQ"^ 



with y being the inelasticity variable. In the laboratory frame, the dominant mechanism is the boson-gluon fusion 
7*G — !■ cc. Hence, the charm structure function is directly driven by the gluon distribution and provides constraints 
for the gluonic function. In leading order (LO), it is written as ]T^, 

Stt Ja,x y ^^"^ \ y 



F2{x,Q\mi)^^^ (-, -f PrG(y,A^^), (28) 



where Oc = (1 -|- -^)- The mass factorization scale lies in the range < < 4 (Q'^ + 4 m^). Such a scale introduces 
an uncertainty of about 10 %, and an additional source of uncertainty is the charm mass, in general ranging on 
1.2 < rric < 1.7 GeV. The standard QCD coefficient function is labeled by Cg 2 {z, 

Experimentally, the latest measurements of the charm structure function are obtained by measuring mesons D* ^ 
production | |3^ . From the theoretical input, it was used NLO coefficient functions, considering charm mass rric = 1.4 
GeV and factorization-normalization scale — + 4 m^. The function F|'^(a:,(5^) shows an increase with 

decreasing x at constant values of , whereas the rise becomes more intense at higher virtualities. The data are 
consistent with the NLO DGLAP calculations. Concerning the ratio R^'^ = F2^ /F2^ the charm contribution to F2 
grow s steeply as x diminishes. It contributes less than 10% at low and reaches to about 30 % for > 120 GeV^ 

Once more the color dipole picture will provide a quite simple description for the charm structure function in a 
factorized way. Now, the Glauber-Mueller dipole cross section is weighted by the photon wavcfunction constituted by 
a cc pair with mass nic- Our expression for the charmed contribution in deep inelastic is thus written as 

Ff{x,Q') / d'r dz {\^¥iz,r)\' + \nHz,r)\') ag'^L{x,r') (29) 

iTT^acm Jo Jo 

where |^'^'^^ (z,r)p is the probability to find in the photon the cc color dipole with the charmed quark carrying 
fraction z of the photon's light-cone momentum with T, L polarizations. For the correspondent wavefunctions, the 
quark mass in Eqs. (^|J|) should be substituted by the charm quark mass rric. Here, we should take care of the 
connection between the Regge parameter x — [W'^ +Q^)/(Q^ + 4:171^) and the Bjorken variable XBj- For calculations 
with the light quarks these variables are equivalent, however for heavier quarks the correct relation is ||29|] , 



13 



C 0.4 
0.2 
0.0 

0.4 
0.2 



0.0 
10 



Q'=1.8 


Q^=4 




\ Q'=18 

1 1 1 . ^TTm, 


\ Q'=30 

1 1 1 


1 1 


, Q'=60 



10 10" 



10" 



10 10" 



10"^ 

X 



10" 



10 10 10 10" 



FIG. 7. The Glauber-Mueller result for the ff*^ structure function as a function of Bjorken variable x at fixed virtualities 
(in GeV^). One uses charm mass nic = 1.5 GeV, target size R\ = 5 GeV ^ and frozen gluon distibution at large r. Data from 
ZEUS Collaboration [33] (statistical errors only). 

In Fig. (|^) we show the estimates for the charm structure function as a function of xb} [Eq. (^0|)] for representative 
virtuahties. In our calculations, it was used the charm mass rric = 1.5 GeV, the target size R\ = 5 GeV~^ and 
the frozen gluon distibution at large r. We have verified a small soft contribution as in the Fl case, decreasing as 
the virtuality rises. There is a slight sensitivity to the value of the charm mass, increasing the overall normalization 
as rric diminishes. Such a feature suggests that the charm mass is a hard scale suppressing the non-perturbative 
contribution to the corresponding cross section, which is in agreement with the recent BFKL color dipole calculations 
of Nikolaev-ZoUer ||2^ and Donnachie-Dosch . 

Regarding the Breit system description, in Ref. |ll3|| it was found strong corrections to the charm structure 
function, which are larger than those of the F2 ones making use of expression (^8|). Considering the ratio 
i?2 = F2'^^{x,Q'^)/F2^'^^^^{x,Q'^), the corrections predicted by the Glauber-Mueller approach would reach to 
62 % at values of ln(l/a;) « 15 (THERA region). Then, an important result is a large deviation of the standard 
DGLAP expectations at small x for the ratio R'^'^ = F^"^ jF^ due to the saturation phenomena (unitarization). With 
our calculation one verifies that it is obtained a good description of data in both reference systems, suggesting a 
consistent estimation of the unitarity effects for this quantity. 



V. CONCLUSIONS 



We study the dipole picture for the description of deep inelastic scattering, focusing on observables driven directly 
by the gluon distribution. Starting from the dipole cross section provided by the Glauber-Mueller approach in QCD, 
we perform estimates for the inclusive structure function F^-, the longitudinal function Fj^ and the charm structure 
function on the proton i^2^^. 

For each of the observables discussed, we obtain theoretical estimates, in the rest frame, without further fitting 
procedure, in good agreement with the updated data from HERA. The resulting calculations corroborate a quite 
consistent picture for the unitarity corrections from the Glauber-Mueller approach in both Breit and rest reference 
systems. In the laboratory frame the unitarity effects are connected with the gluon distribution function, whereas in 
the color dipole framework the basic block is the dipole cross section which is corrected considering saturation effects. 

The small transverse separation r region is dominated by the leading log DGLAP formalism, with the additional 
ingredient of unitarization phenomenon as the momentum fraction acquires quite small values. Such corrections 
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are associated with the taming of the gluon distribution in the very small x region, in general named saturation 
regime. However, it should be stressed that the Glauber-Mueller approach and similar eikonal-like models provide 
a logarithmic ln(l/a;) asymptotic behavior for the inclusive structure function and gluon distribution, instead of a 
constant value for asymptotic energies. 

The large transverse separation is described by non-perturbative aspects of QCD. Since this domain is not well 
determined at the moment, some modelling of the soft region is needed. In this work we choose the ansatz in which the 
gluon distribution is frozen for virtualities above a cut radius > r^^^ , which corresponds to the region < Qq. A 
convenient choice for the gluon pdf in order to cover the widest possible kinematical window diminishes the uncertainty 
coming from the soft sector. The most appropriated input is the GRV94 parametrization, where Tcut = 0.6 fm is 
found, whereas it can take values Tcut = 0.4 — 0.5 fm for the more recent pdf's. Throughout the paper we used the 
target size R\ — 5 GeV^^, which corresponds to strong unitarity corrections. 

When performing a comparison with the phcnomcnological model GBW, we have found that the Glauber-Mueller 
approach underestimate the dipole cross section from GBW at not small x > 10~^. Instead, for very low x < 10~^ 
the Glauber-Mueller overestimates GBW due to the strong increasing of the gluon function in this region. Concerning 
the saturation model with DGLAP evolution (BGK), it is a particular case of the GM approach when considering 
central scattering 6 = 0. Despite that the BGK model matches DGLAP evolution, Glauber-Mueller describes more 
properly the realistic impact parameter dependence of the process. Moreover, in GM the extension to the nuclear 
case is built in. 

When considering the structure function F2, we have found that it is dominated by small transverse distances 
contributions. However, a non-negligible content from the soft sector is present. Moreover, the photon wavefunctions 
enhances the dipole cross section into smaller dipole sizes, since the weight function selects smaller r as the virtuality 
Q- diminishes. Our estimates here are parameter-free, however a fine tunning of the parameters can improve the 
data description. Furthermore, we notice that in calculations from |10|, only the aligned jet dipole configuration z, 
(1 — z) « (and only transverse contribution) is considered, whereas we take into account all configurations, including 
the symmetric ones. Thus, all dipole sizes, even those from the non-perturbative region are included in our results. 

Concerning Fl, the estimates are consistent with the previous calculations in the Breit system and are in good 
agreement with data. A remarkable feature is that the Glauber-Mueller approach in the color dipole framework gives 
important higher twist contributions to the leading twist calculation in a simple way. As is well known, is the 
main quantity to study the expected higher twist effects in low virtualities. 

The function is directly dependent of the gluon distribution and important unitarity corrections had been found 
when considering the Breit frame. Here, we verify consistent results in the rest frame in comparison with the previous 
ones in the fast proton system. We verified that the charm mass suppressed soft contributions in comparison with 
the F2 case, and the results present a slight dependence with the specific value of TOc- 

In conclusion, the Glauber-Mueller approach provides a well stablished formalism to take into account the unitarity 
effects. It allows to estimate the higher twist contributions to relavant observables in a simplified way, while matching 
DGLAP evolution equation at Born level and including the impact parameter dependence properly. 
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